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Density functional theory (DFT) is used to study vibrations, electrical dipole moments, and 
polarizabilities of NaF clusters. Because of prior experimental and theoretical studies, this is a good 
model system for tracking the evolution of the properties from diatomic molecule to bulk crystal. 
The ratio of vibrational to electronic contributions to the polarizability increases dramatically with 
size N in the closed shell clusters NajvFjv. The open shell system Nai4Fi3 has a greatly enhanced 
electronic polarizability. Contrary to previous studies on this system which treated only the outer 
electron by quantum mechanics, we find the Oh cubic structure to be stable relative to the polar 
distorted structures such as C3,,. 



I. INTRODUCTION 

Rayane et al. 1 have pioneered the experimental study 
of dipole moments and polarizabilities of NaF clusters, 
focussing on systems such as NajvFjv-i with one electron 
outside a closed shell. However, even the closed shell 
systems Na^vF^ have dielectric properties of considerable 
interest. A study of this can illuminate such issues as how 
ferroelectric dipolar order of the bulk evolves up from 
nanocrystals. Our main focus is the quantity ajv, defined 
as the polarizability per molecular unit of a cluster of N 
molecular units. 



a N ,a/3 = -^dfia/dFp 



(1) 



where jl is the total dipole moment of a cluster or crystal 
and F is the external electric field. 
The Clausius-Mossotti relation 



cvoo = lim a N = — 

N^oo 1 



Oil 



(47rn/3)ai 



(2) 



is often used to relate the crystal polarizability per 
molecule a^, = (e— l)/47m to the molecular polarizabil- 
ity a\ and the molecular density n. This equation is ex- 
act for a bulk crystalline solid where point-polarizable 
and isotropic molecules sit at sites of cubic symmetry. 
It takes into account the local electric fields of all other 
induced dipoles which enhance the applied field at the 
cubic sites. This provides the simplest model for a fer- 
roelectric instability when ot\ is large enough. Numeri- 
cal linear algebra allows the formula to be extended to 
less symmetric systems of point-polarizable molecules 2 . 
However, no solid sequesters all polarizable entities accu- 
rately into point objects, or has only dipolar interactions, 
so the applicability of the formula to real ionic crystals, 
while traditional, can be regarded as largely wrong 3 , or 
alternately as basically right but only after the ions are 
altered by crystal field and confinement effects 4 . Here 
we use electronic structure theory (density functional 



theory, or DFT) to provide "data" on the evolution of 
polarizability 5 . 

Semiconductors like Si 6 and CdSe 7 have received a 
large amount of attention. The alkali halides, however, 
provide simpler model systems. Hudgins et ai. 8 have doc- 
umented experimentally the stability of simple cuboidal 
fragment structures in NaCl. Much modeling has been 
done, mostly using classical models such as the shell 
model 9-14 . Electronic structure methods have been ap- 
plied to diatomic alkali halide molecules 15-17 and various 
alkali halide clusters 18-20 . 

We choose NaF as a model system for studying the evo- 
lution toward bulk properties as well as the idiosyncrasies 
of individual nanoparticles. The bulk 21 and surface 22 vi- 
brational properties are both known. Foremost among 
idiosyncrasies of nanophases is the fact that non-zero 
equilibrium electrical dipoles fx are necessary, just for 
reasons of cluster geometry and symmetry. This is well 
known 23,24 in materials like CdSe whose bulk phase lacks 
inversion symmetry. It is less well known that dipoles 
must occur even in rather symmetrical cuboid rocksalt 
fragments, unless they have high enough symmetry (T^ 
or symmetries which include inversion) to forbid a dipole. 

Concerning electrical polarity of nanoclusters, there is 
a large literature of experimental and theoretical stud- 
ies for metal clusters, especially simple metals 25,26 (e.g. 
Na). There are also studies of polarity in semicon- 
ducting nanoclusters 27 . Microscopic dielectric screen- 
ing theory has been developed and implemented for 
semiconductors 28,29 . In ionic nanophase materials no 
systematic study of polarity using full electronic struc- 
ture theory is known to us. One new feature in ionic 
materials, not seen in metals or homopolar semiconduc- 
tors, is a large contribution to polarity coming from nu- 
clear displacements. We call this contribution the "vi- 
brational" part a v ib- Heteronuclear materials with par- 
tially ionic nature, such as the rocksalt crystal GeTe, 
can have ferroelectric instabilities driven partly by vi- 
brational polarization 30 . 

NaF is a good choice for a first study for two rea- 
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TABLE I: Experimental and theoretical properties of the NaF 
diatomic molecule. 



present calc. 


expt. expt. 


theory 


theory 


NRLMOL 


ref.33 ref.34 


ref.17,19,20 


ref.15,16 



d (ttfi) 
[i (ea B ) 
uj (cm -1 ) 
a e i,± (A 3 ) 
aei,|| (A 3 ) 
avib,|| (A 3 ) 



3.69 
3.07 
510 
2.39 
3.33 
0.86 



3.64 
3.21 
536 



3.64 



536 



3.6 
3.22 
540 
1.82 
1.95 
1.04 



3.2 



TABLE II: Experimental and theoretical polarizabilities (in 
A 3 ) of Na and F atoms and ions. SCF means self-consistent 
Hartree-Fock, and ACPF (average coupled pair functional) 
includes correlation. 





Na 


Na+ 


F 


F" 


NRLMOL 


23.18 


0.144 


0.536 


1.23 


expt. ref 35 


24.11 








SCF ref. 16 




0.139 




1.58 


ACPF ref. 16 




0.148 




2.84 



sons. First, there is a lot of previous work on struc- 
tural properties of this and related alkali halide clus- 
ters, for which structures are generally simple fragments 
(nanocrystals, or "cuboids") of the bulk rocksalt struc- 
ture. Second, because of the small polarizability of the 
F~ ion, and the similar sizes of Na + and F~ ions, NaF 
should be simpler than all other alkali halidcs, potentially 
more amenable to analysis using models in the spirit 
of Clausius-Mossotti. There have been surprisingly few 
calculations using full electronic structure theory; non- 
stoichiometric varieties (which we largely ignore) have 
received the most attention 31 . We postpone the study 
of classical pictures for a future paper, and present here 
the results of a quantum study using DFT. This paper 
is based largely on the results of the MA thesis of C. 
Schmidt 32 . 



II. DIATOMIC MOLECULE 

The NaF diatomic molecule or "monomer" allows a 
good comparison between our calculations and both 



experiment 



33-35 



and previous theory 



16,19,20,36 



We use 



the package NRLMOL 37 which allows versatile and ac- 
curate DFT calculations for molecules. We choose the 
PBE 38 version of generalized gradient approximation 
(GGA), and include a spin density functional for systems 
with uncompensated spins {i.e. the Na and F atoms and 
Nai4Fi 3 .) Results are shown in Table I. We also show 
in Fig. 1 the total energy from DFT relative to isolated 
atoms. The experimental dissociation energy of 5.3 eV 39 
is quite accurately reproduced, as are the bond-length, 
dipolc moment, and vibrational frequency. 

Can the diatomic results be related to atomic, ionic, 



or bulk properties? Computed atomic and ionic polariz- 
abilities are shown in Table II. The sum of DFT polariz- 
abilities of the Na + and F _ ions is 1.37 A 3 . This is close 
to the polarizability cvoo ~ 1.45 per formula unit of the 
solid. However, Clausius-Mossotti theory gives the value 
ai=1.16 A 3 from the crystal data. It is difficult to know 
how to interpret these numbers. Fowler and Madden 4 
make three relevant observations. (1) Correlation correc- 
tions greatly enhance the polarizability of the larger F - 
ion, doing little for the isoclcctronic but smaller Na + ion. 
(2) Putting the ion in a crystalline environment makes 
the F~ ion shrink, both because of the attractive crystal 
field of the surrounding ions and because of the repulsion 
of the surrounding electrons. The resulting smaller ion 
has a polarizability which is much less altered by corre- 
lation. (3) The resulting smaller polarizability estimated 
for the "ion-in-a-crystal" is perhaps quite consistent with 
Clausius-Mossotti. It is unfortunate that there are no di- 
rect experimental data for the polarizatility of the ions 
Na + or F~ . A long list of highly diverse values for the F~ 
polarizability is tabulated by Iwadate and Fukushima 40 . 
None are direct measurements. 

It is difficult to know how to relate these ideas and 
numbers to our calculated polarizability of the diatomic. 
In zeroth approximation the equilibrium dipole of the 
diatomic is ed and the electronic polarizability is the 
sum of the ionic polarizabilities, or 1.37A 3 . The com- 
puted dipole is smaller by 17%, but the computed po- 
larizability is far larger than the zeroth order estimate. 
Using the average a = (a xx + a yy + a zz )/3, the diatomic 
has <5i iC i=2.70, larger by 2 than the sum of free ions. 
Also the diatomic has anisotropy Aa/a=34%, where 
Aa 2 = [(a xx - a yy ) 2 + (a yy - a zz ) 2 + (a zz ~ a xx ) 2 }/2. 
In an attempt to account for these discrepancies, three 
corrections were tested. (1.) Each ion polarizes in the 
self-consistent field of the other, and these local field 
corrections are accounted by solving two coupled linear 
equations for the moments fii and \i<i. Unlike Fowler and 
Madden 4 , we did not correct the polarizability of the 
ions for their altered environments. The corrections de- 
pend on the ratios <j)± = a±/d 3 for the two ionic species, 
where d is the internuclear separation. For Na + <p + = 
0.0193, and for F" </>_ = 0.1652. The reduced ratio 
0rcd = 4>+4>-/( ( f ) + + 4>-) is 0.0173. The corrected val- 
ues of the permanent dipole 41 and polarizability tensor 15 
are 



1 - 0+ - 0_ 



/! = ed 
a\\ = (a + + a_) 
a± = (a + +a_) 



1 - 40+0_ 
" l + 4c/> red 



- 2<^red 



= 0.80ed 
1.08(a+ + a_) 
= 0.97(a+ + 



(3) 

The result of simple screening of the naive dipole ed by 
point ionic polarizability agrees closely with the com- 
puted total dipole (Table I), but the enhancement of a 
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1 2 3 4 5 
separation (Angstroms) 



ity, we look instead to the bulk solid for guidance, and 
here a surprise occurs. The static etot = £dc and optical 
£ei = £ opt dielectric constants of NaF are 5.1 and 1.74 ; 
the effective polarizabilities a = (e — l)/47m are 8.04 
A 3 (electronic plus vibrational) and 1.45 A 3 (electronic 
only). Thus the vibrational contribution is 5 times the 
electronic contribution to the bulk polarizability, whereas 
in the diatomic molecule according to Table I the ratio 
is 0.29/2.70, 50 times smaller (this uses the average of 
the diagonal elements of the a tensor.) An alternate 
analysis of the bulk is to extract "molecular" polariz- 
abilities from the Clausius-Mossotti formula. These are 
1.16 A 3 (electronic) and 2.29 A 3 (vibrational). The ra- 
tio is reduced from 5 to 2 but is still 20 times larger 
than in the diatomic molecule. This remarkable differ- 
ence comes partly because the electronic polarizability is 
reduced (per molecule) in the bulk, but mostly because 
the vibrational part is increased. 



FIG. 1: Density functional total energy of NaF diatomic 
molecule relative to free atoms. 

parallel to the molecular axis, and smaller enhancement 
perpendicular, while in the right direction, are much too 
small. It appears that when higher-order multipoles are 
included, they make corrections comparable to the dipo- 
lar ones so the agreement with the total dipole mo- 
ment is partly accidental. (2.) Nonlinear corrections 
are not negligible since the local fields are large. Since 
the ions are inversion-symmetric, the leading correction 
is [i = (a + f3F 2 /2)F. Using NRLMOL we computed 
(3—7.9 and 0.06 for F~ and Na + , in atomic units (a^/e 2 .) 
The ions sit in fields F = e/r 2 which are large, and the 
non-linear correction is estimated to be +12% for the 
polarizability of the F part of NaF, and +1% for the Na 
part; again the direction is right but the effect is much 
too small. (3.) There is incomplete charge transfer be- 
tween Na and F in the diatomic. Using the computed 
zero-field dipole as a measure, the static charge \ijd is 
±0.83e, which suggests that the Na could be treated as 
83% Na + and 17% Na. Because of the giant polarizabil- 
ity of neutral Na (see Table II), it is easy to rationalize 
a large correction, but less easy to make a simple clas- 
sical model for both tensor components of the electronic 
polarizability. 

The vibrational part of the polarizability is a new effect 
that occurs in the diatomic and has no counterpart in 
atoms or ions. We compute it from the formula 43 ' 44 

a v ib, a p = ^ Z a .i^(K~ 1 ) i ^.j lJ Z'^j 3 (4) 

where Z a ^^ is the dynamic effective charge tensor 
d^a/dui^, is the /i-th Cartesian component of the 
displacement of the i-th atom, and K is the force con- 
stant tensor. 

Since atoms and ions have no vibrational polarizabil- 



III. CLUSTERS 

We have studied Na^F^ for N=2, 4, 6, 9, 18, and 32. 
For contrast, we also studied the system Nai4Fi 3 which 
is a symmetric 3x3x3 cuboid of Oh symmetry. Except 
for N—2, and the ring structure of N=4, these are all 
cuboids, shown in Fig. 2. Vibrational frequencies and 
infrared activities used to compute a v ib are summarized 
in Fig. 3. Polarizabilities are summarized in Figs. 4 and 
5, and polarizabilities, formation energies, and perma- 
nent dipole moments are summarized in Table III. 



A. Na 2 F 2 

This cluster is not a cuboid, but is included because 
it has previously been studied experimentally 46-49 and 
theoretically 9 ' 19,20 . Our DFT gives a binding energy rel- 
ative to two separated diatomics of 2.4eV. Previous cal- 
culations by Dickey et al. 19 gave 2.6cV, which agrees well 
with the experiment of Eisenstadt 46 . The geometry is a 
rhombus with Na-F-Na angle 94.7° and Na-F bond length 
3.96 as, agreeing with Dickey et al. and with experi- 
ments by Hartley and Fink 47 and Lapshina et al. 4& . 

The DFT prediction of polarizabilities is given in Table 
III. The electronic contribution is surprisingly isotropic. 
Per ion pair, the polarizability a c i.2, plotted on Fig. 
4 is significantly reduced compared with the diatomic 
molecule, and closer to the bulk limit. The vibrational 
polarizability has a large anisotropy dominated by out 
of plane motions, and per ion pair, a v ib,2 is much larger 
than the diatomic and out-of-line with the results for 
larger clusters, already close to the bulk value. The 6 
vibrational normal modes consist of three even, Raman 
active modes (A g at 205 cm -1 , Bi g at 326, and another 
A g at 377) and three odd, infrared active modes, roughly 
equally intense (Bi„ at 149 cm -1 , B^ u at 362, and B2« 
at 373). In assigning these symmetries, we used a co- 
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FIG. 2: Structures of NaF rocksalt fragments (cuboids, ex- 
cept for Na2F2) given by DFT. The Na4F4 ring structure 
and Oh 4x4x4 cuboid Na32F32 are not shown. The 
NagFg and NaisFis C 4 i, structures have even numbers of 3x3 
planes; planes are charged alternately + or -, so the fragments 
are polar. The Nai4Fi3 Oh cuboid has an odd number of 
3x3 planes and is non-polar. It consists of a closed shell 
positively charged cluster plus a compensating electron dis- 
tributed mainly on corner Na + ions. 



ordinate system where the x-axis passes through the Na 
atoms and the y axis through the F atoms. These num- 
bers can be read from Fig. 3, which also shows the rela- 
tive infrared intensities of the active modes. Our frequen- 
cies agree to better than 10% with unrestricted Hartree- 
Fock calculations 19 and somewhat less well with shell 
model calculations 9 . A matrix- isolation experiment 49 
has reported two Na2F 2 infrared active vibrations above 
a spectrometer cutoff at 190 cm -1 . The experimental 
numbers 360 and 363 cm" 1 are close to our values. It 
is the low frequency B lu mode, polarized out of plane, 
which is responsible, because of its low frequency, for the 
large zz-component of vibrational polarizability in our 
calculation. 



ation and if harmonic approximation is sufficient, and if 
we ignore rotations, then statistical mechanics says that 
the fraction of rings fn should be given by 



l/f R = 1 + [] 



sinh(/37iu;f/2) 
smh{l3hujf /2) 



(5) 



where Lof is the z-th vibrational mode of the ring and 
is the i-th vibrational mode of the cuboid. At low 
T = l/fikB, the ring population is exponentially small, 
but as T increases, the lower vibrational frequencies give 
an entropic stabilization of the rings. In the classical 
limit pUco/2 <C 1, the temperature Tr at which the ring 
population is 1/2 is 



fcsT^A/^logtaf/wf) 



(6) 



Our computed vibrational spectra, shown in Fig. 3, give 
the logarithmic denominator to be 8.5, greatly reducing 
the temperature of isomerization. This reduction fac- 
tor should be slightly increased to take into account the 
higher rotational moments of inertia, and thus higher ro- 
tational entropy, of the rings. 

The polarizabilitics, shown in Table III, are surpris- 
ingly isotropic for the ring. Electronic parts are surpris- 
ingly similar for the ring and cube geometries, while the 
vibrational polarizability is bigger by a factor 4 in the 
ring because of the low frequency infrared activity shown 
in Fig. 3. 

An interesting property of the Na 4 F 4 T^ structure is 
that the ground state has an octupole (xyz) moment, but 
no quadrupole moment. In the point group T^, a vector 
(x,y, z) and an off-diagonal second-rank symmetric ten- 
sor (yz, zx, xy) have the same symmetry. This means 
that an applied electric field F z z induces both a dipolc 
moment \i z and a quadrupole moment ©^j, 50 . We calcu- 
lated the octupole moment, which is equivalent to static 
point charges of ±0.866e on the ions. This can be com- 
pared with the static charge ±0.83e deduced from the 
dipole of the monomer NaF. The induced quadrupole was 
computed to be 



Q xv = lJ d 3 rxyp(r) = (-8.2a%)F z . 



(7) 



B. Na 4 F 4 

This cluster has a relatively low-lying structural iso- 
mer with a D^h rmg geometry which has been studied 
theoretically 14 . We find the Td cuboid geometry to be 
lower in energy than the ring by A=0.75 eV, twice as 
much as the figure 0.37 eV mentioned by Calvo 14 . The 
vibrations of the ring have a wider range of frequencies 
than any of the cuboids, with many low- frequency modes. 
We estimate the temperature of isomerization to be 1000 
K in the following way. If only two isomers need consider- 



C. Larger Stoichiometric Cuboids 

The vibrational eigenfrequencies are shown in Fig. 3. 
The polarizabilities are shown in Fig. 4, and contrasted 
with crystalline compounds in Fig. 5. We see a non- 
monotonic but fairly smooth evolution of both vibra- 
tional and electronic polarizabilities with linear dimen- 
sion of the cluster, provided we ignore the planar species 
Na 2 F 2 and D 4h Na 4 F 4 . The Na 4 F 4 and Na 34 F 32 T d 
structures are the first two in the series of T^ structures 
of 2n x 2n x 2n atoms with isotropic polarizabilities whose 
properties should scale smoothly and perhaps rapidly to 
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TABLE III: Polarizabilities (in A 3 ), formation energy E{ (in eV), and permanent dipole moment y, (in atomic units ea B ) for 
NaF molecules and nanocrystals. 
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49.4 


49.4 


89.1 


89.1 


89.1 
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"Ref. 39 
'Ref. 46 



the bulk limit. The electronic polarizability does ap- 
proach the bulk limit rapidly, but the vibrational part 
does not. There is a subtlety in the vibrational polar- 
izability of the bulk. The Lyddane-Sachs- Teller relation 
shows that the splitting of TO and LO optical phonon fre- 
quencies is tied to the splitting e^c — e pt- In point group 
Td or Oh, the infrared optic modes belong to triply de- 
generate representations - there is no distinction between 
LO and TO. The bulk crystal violates this symmetry only 
because the distance of spatial variation of the probing 
electric field is smaller than the sample size. Therefore 
the bulk limit for a v ib may not be reached until the sam- 
ple size exceeds the wavelength of the radiation used to 
measure £dc- 

We have done three clusters with uniaxial anisotropy. 
NagF 6 is first in a sequence of long skinny cuboids of 
2 x 2 x (2n + 1) structures with I)2h symmetry or of 
2 x 2 x (2n + 2) structures with T)2d symmetry. Actually 
the former group is not strictly uniaxial, but has weak 
splitting of the equatorial two-fold axes, which accounts 
for a very weak splitting of the xx and yy parts of the a 
tensors. NagFg and Na^Fis are the first two of the 
sequence of 3 x 3 x 2n structures. The polarizability is 
relatively isotropic, except that the vibrational polariz- 
ability of NaigFis is 20% enhanced along its long axis. 
The clusters with dimensions 3 x 3 x (2n + 1) have an 
odd atom count and are thus non-stoichiometric. The Oh 
3x3x3 cluster Nai4Fi 3 is discussed below; the higher 
clusters have D 4 /j symmetry. 



D. Permanent Dipole Moments 

The two cuboids of symmetry, NagFg and 

NaisFig, have permanent dipole moments just of struc- 
tural necessity. Naively, if the cuboid were an exact 
sodium chloride structure fragment, with the NaF bulk 
lattice constant of 4.62 A, and the charge of ±e were as- 
signed to the ions, then the dipole moments would be 2.31 
eA and 4.62 eA for the smaller and larger cuboids. How- 



ever, the ions have somewhat distorted positions, and the 
structural dipoles are reduced to 1.34 eA and 3.16 eA re- 
spectively, assuming unit quantized charges. The actual 
DFT values, shown in Table III arc 1.06 eA and 2.71 eA, 
respectively. These numbers show that the total dipole is 
reduced by screening, both "vibrational" and electronic. 
Vibrational screening is larger than electronic, in line 
with the fact that vibrational polarizabilities are roughly 
twice electronic polarizabilities. The slab-shaped NagFg 
is more affected by screening than the needle-shaped 
NaisFis, as one expects from the well-known depolar- 
ization effects in slab vs. needle shaped samples. The 
dipole of Nai 8 Fi 8 was previously calculated by Rayane 
et al. 1 to be 7.5 D = 1.84 eA, significantly smaller than 
our value. Their calculation was part of an investigation 
of non-stoichiometric Na^F^r-i species such as NaigF^. 
Only the outer electron was treated by quantum mechan- 
ics. The stoichiometric species has no outer electron, so 
apparently their number is based on a classical polariz- 
able ion model, which apparently overscreens. 



E. Non-stoichiometric O h Nai 4 Fi3 

This is such an interesting cuboid that we could not 
resist departing from our main theme of polarizability 
evolution in stoichiometric cuboids. The last electron re- 
sides outside the closed shell of Na^FT^. Alkali halidcs 
with one electron outside a closed shell have attracted 
attention 1,51 since the prediction by Landman et al. that 
such clusters might have a second-order Jahn- Teller in- 
stability. This, and the subsequent studies by Durand et 
a l 53,54 kgpt only a single electron in the quantum part 
of the dynamics, while Ochsenfeld and Ahlrichs 55 have 
studied several alkali halides using all-electron theory. 

The single outer electron gives rise to a giant electronic 
polarizability shown in Table III. On a per molecule ba- 
sis, the total polarizability of 313 A 3 , shared by 13.5 
molecules, is larger by 15 than the interpolated value 
~ 1.5 seen in Fig. 4. On an absolute basis, it is larger by 




o electronic 
□ vibrational 



100 200 300 400 500 
normal mode frequency (cm ) 



FIG. 3: Integrated vibrational densities of states for various 
clusters of NaF. The degeneracy of each mode is coded in the 
step increase, as the spectrum steps upward to the final mode 
3N — 6. Vertical bars are the relative infrared activities of the 
various modes; the absolute scale is not given. 



13.5 than the polarizability of the neutral Na atom. On a 
per molecule basis, the vibrational polarizability of 24.7 
is smaller by 0.6 than the interpolated value of the sto- 
ichiometric clusters. Apparently the loose last electron 
is exceedingly polarizable, and partly screens the change 
in dipole per unit displacement that causes vibrational 
polarizability according to Eq. (4). 

The large polarizability gives one a right to expect that 
Nai4F 13 might distort to a lower symmetry polar struc- 
ture such as C3t,, as was discovered theoretically in a 




(1/N) " N=number of molecules in cluster 



FIG. 4: Scaling with l/iV 1/3 of the polarizability per NaF 
molecule a at of clusters of various sizes NajvFjv. The di- 
atomic or monomer (N = 1) and the dimer (N = 2) are not 
fragments of rocksalt, but the other clusters, iV=4,6,9, and 
18, are fragments of the bulk crystal (N = oo). The straight 
lines just connect the diatomic to the bulk as if the corrections 
were linear in 1/N 1 ^ 3 . 



restricted model by Landman et al. 52 in the Nai4Cli3 
cuboid. In a naive view, point-polarizable ions have a 
density-dependent critical a i,c beyond which ions po- 
larize spontaneously. In cubic crystals it is 3 /Aim which 
equals 5.88 A 3 for the density n in bulk NaF. In fragments 
a c approaches the bulk limit rapidly 2 . The per molecule 
electronic polarizability of Nai4Fi3 is way beyond that 
limit, and thus unstable in this model. However, the gi- 
ant polarizability is almost all associated with the delo- 
calized last electron, making the Clausius-Mossotti anal- 
ysis inappropriate. 

The highest occupied molecular orbital (HOMO) is 
non-degenerate but lies not far below other levels with 
which it is forbidden by symmetry to mix. The second- 
order Jahn- Teller effect is a distortion motivated by 
the "desire" for such a mixture in order to lower 
(quadratically in distortion amplitude) the HOMO en- 
ergy. Kristoffel and Konsin 56 argue that such a mech- 
anism is at work in the ferroelectric distortion of GeTe 
crystals, involving the mixing between states just below 
and just above the narrow gap of this IV- VI semiconduc- 
tor with two more electrons than are needed to fill the 
usual s — p shell. The instability of Landman et al. 52 is 
a nanocrystalline analog. 

We find that this instability does not happen in 
Nai4Fi3. This contradicts other theories and the inter- 
pretation by Rayane et al. of their experiment. Theories 
which keep only the extra outer electron are not reliable 
enough if the energy balance is delicate. Rayane et al. 1 
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FIG. 5: Log-log plot of vibrational polarizability (vertical) 
vs. electronic polarizability (horizontal). Open squares are 
the NaF cuboid clusters. Closed circles are from crystal data. 



found a very large deflection of neutral Nai4Fi3 parti- 
cles in electric field gradients. There was a moderate 
T-dependence, but not the full 1/T of the usual Debye- 
Langevin polarizability /i 2 /3fcsT of a molecule with a 
permanent dipolc \x. Nevertheless, the authors believe 
their data to be consistent with a permanent dipole, al- 
though they do not claim that their analysis is firm. Our 
finding is that the electronic polarizability is indeed ex- 
tremely large and will give large deflections, but there 
is no symmetry-breaking or permanent dipole moment. 
We believe the data may be consistent with this alternate 
interpretation. 



IV. DISCUSSION 

A. Formation energy 

We define the formation energy Ef as the energy per 
ion pair relative to the energy of separated neutral atoms. 
Values of Ef are given in Table III. The increasing for- 
mation energy with cluster size can be organized phc- 
nomenologically according to the numbers of ions with 
various geometric sites. We distinguish corner, edge, 
face, and interior ions. In the Na^F^ cuboids these sites 
occur as pairs of oppositely charged ions. For example, 
Na4F 4 has 4 corner pairs, Na 6 F 6 has 4 corners and 2 
edges, NagF 9 has 4 corners, 4 edges, and one face, and 
NaigFig has 4 corners, 8 edges, 5 faces, and 1 interior ion 
pair. The values of Ef then correspond to the assignment 
7.01 eV for corner pairs, 7.53 eV for edge pairs, 7.63 eV 
for face pairs, and 9.11 eV for interior pairs. This last 
is close to the experimental bulk value of 9.30 eV. These 



assignments then predict Ef for larger cuboids. Na32F 32 
has 4 corners, 12 edges, 12 faces, and 4 interior pairs, 
so is predicted to have Ef— 7.70 eV, very close to the 
computed DFT value of 7.6 eV. 



B. Trends of Polarizability 

The trend to smaller electronic polarizability as size 
increases should certainly be attributed to enhanced sta- 
bilization of the filled p shell as the ion acquires larger 
numbers of neighboring counterions. The HOMO and 
LUMO states of a finite cluster are located mostly at the 
surface. The HOMO-LUMO gap is reduced in clusters 
relative to bulk, corresponding to a cluster analog of sur- 
face localized states on a bulk surface. As the cluster 
size increases, electronic spectral weight gets pushed up- 
ward toward the bulk band gap of lleV 58 . In bulk alkali 
halides, an excess electron is a delocalized carrier, but an 
excess hole self-traps on a dimerized anion pair (V^ cen- 
ter). On clusters, one would expect the trapped hole to 
lie on surface anions. The excess electron, after migrat- 
ing to the surface, might also self-trap. This would be a 
"pseudo-Jahn- Teller polaron" . The fact that the second- 
order Jahn- Teller effect does not occur in NaigF^ (ac- 
cording to us) suggests that the surface electron might 
not trap but remain delocalized. We note that localiza- 
tion of the excess electron could be enhanced within a for- 
malism where the electronic self-interaction (SI) is fully 
corrected 59 . Further investigation of this effect within 
an Si-free framework would be necessary to rule out the 
possibility of surface trapping. 

The enhanced vibrational polarization of larger clus- 
ters is harder to explain. Large electronic polarizabil- 
ity should diminish vibrational polarizability by screeing. 
The dynamical effective charge matrix (change of electri- 
cal dipole per unit displacement) is reduced by the po- 
larizability of the ions which are displacing (see Eq. 3.) 
We find little support for this idea in the trends of di- 
electric constants of other alkali halides, shown in Fig. 5. 
Here the value of a per molecule is (e — l)/47rn, with the 
optical dielectric constant giving the a e i as before. The 
vibrational polarizability seems roughly independent of 
the electronic polarizability. 
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